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ABSTRACT 

We study the hydrostatic density structure of the inner disc rim around Herbig Ae stars using 
the thermo-chemical hydrostatic code ProDiMo. We compare the Spectral Energy Distri- 
butions (SEDs) and images from our hydrostatic disc models to that from prescribed density 
structure discs. The 2D continuum radiative transfer in ProDiMo includes isotropic scat- 
tering. The dust temperature is set by the condition of radiative equilibrium. In the thermal- 
decoupled case the gas temperature is governed by the balance between various heating and 
cooling processes. The gas and dust interact thermally via photoelectrons, radiatively, and via 
gas accommodation on grain surfaces. As a result, the gas is much hotter than in the thermo- 
coupled case, where the gas and dust temperatures are equal, reaching a few thousands K in 
the upper disc layers and making the inner rim higher. A physically motivated density drop 
at the inner radius ("soft-edge") results in rounded inner rims, which appear ring-like in near- 
infrared images. The combination of lower gravity pull and hot gas beyond ~ 1 AU results 
in a disc atmosphere that reaches a height over radius ratio z jr of 0.2 while this ratio is 0. 1 
only in the thermo-coupled case. This puffed-up disc atmosphere intercepts larger amount of 
stellar radiation, which translates into enhanced continuum emission in the 3-30 /im wave- 
length region from hotter grains at ~ 500 K. We also consider the effect of disc mass and 
grain size distribution on the SEDs self-consistently feeding those quantities back into the gas 
temperature, chemistry, and hydrostatic equilibrium computation. 
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1 INTRODUCTION 

Circumstellar discs result from the conservation of angular momen- 
tum of a collapsing cloud that forms a star at its centre. Discs of 
various masses and shapes are observed by direct ima ging or spec- 
troscopically from low mass brown dwarfs (e.g.. [Mohantv et al.l 
|2004 l lavawardhana et alj|2003l : iNatta & TestfcOOll) to massive O 
stars iCesaroni et alj 2007b . Spectral Energy Distributions (SED) 
are dominated by dust grain emission and are influenced by disc ge- 
ometry, dust mass, and grain properties. Herbig Ae stars (2—4 Mq) 
are among the most studied pre-main-sequence stars surrounded 
by discs because of their isolated nature and intrinsic brightness 
jNatta et alj|2007l) . Their SEDs can be classified into two groups 
accordi ng to the ratio be tween the far-infrared to near-infrared 
fluxes JMeeus et ai] |200lh . Group I objects show much stronger 
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far-infrared to near-infrared fluxes than their group II counterparts. 
Early versions of the so-c alled puffed-up rim model succeeded 
to ex plain the dichotomy (Dullemond et al. 20ol]; IPominik et"al] 
2003). Group II objects are discs that absorb most of the stellar 
flux within a geometrically thin but optically thick inner rim that 
faces directly the star, depriving the outer disc of radiation. As a 
consequence the cold outer disc is geometrically flatter and emits 
only weakly in the far-infrared. On the contrary, the rim in group 
I objects does not block enough stellar radiation and the outer disc 
flares, emitting strongly in the far-infrared. 2D dust radiative trans- 
fer models that include dust scattering and hydrostatic disc struc- 
ture assuming that the gas temperature is equal to the dust tem- 
perature, which is derived from dust energy b alance, underesti- 
mate the flux in the near-infrared around 3 ftm JMeiier et al ]|2008l : 
IVinkovic et alj2006l) . lvinkovic e t al. (2006) invoke the presence of 
an optically thin halo above the inner disc to explain the excess 
3/im emission. The problem arises from the insufficie nt height of 
the inner rim: a low rim results in a small emitting area. lAcke et al.l 
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(2009) require that their inner rims have to be 2 to 3 times higher 
than the hydrostatic solution that fits the SEDs using a Monte- 
Carlo radiative transfer code. Models that incorporate detailed 
physics of grain evaporation with grain size distribution change 
the structure of the inner rim and help to fit si multaneously the 
SEP and near-infrared interferometric dat a (e.g., Isell a et al]|2006l : 
Tannirkulam et alj2007l ; lKama et al.l20 09). The structure of the in- 
ner rim impacts the overall shape of the SEDs up to a few tens of 
microns. Initially discovered in the SED of Herbig Ae discs, inner 
rim emission is also de tected in discs around classical T Tauri stars 
dMuzerolle et alj2003l) . Observationally, the inner disc geometry is 
varied. Studies with sufficient baseline-co verage show that either an 
extra optical thin c omponent is required I Tanni rkulam et al1l2008l ; 
iBenistv et al.l2010h or that the rim has a skewed asymmetric geom- 
etry dKraus et alj2009h . A detailed modelling of these observations 
is beyond the scope of this paper. 

All previous studies of inner disc rims share the assumption 
that the gas and dust temperatures are equal to simplify the radia- 
tive transfer. In these studies, the temperatures are derived from 
the dust energy balance. This assumption is valid in the inner opti- 
cally thick parts of disc midplanes where gas and dust are thermally 
coupled via in elastic collisio ns, a phenomenon called thermal ac- 
commodation dTielensll2005l) . In all other parts, the gas is mostly 
heated by hot photoelectrons from dust grains and Polycyclic Aro- 
matic Hydrocarbons and attains higher temperatures than the dust. 
The inner rim is naturally higher than when the gas and dust tem- 
perature ar e equal and the overall disc surfaces are more vertically 
extended ( Woit ke. Kamp & Thii r2009 ) . In discs where the gas and 
dust are well mixed (i.e. without dust settling), the change in geom- 
etry can affect the SEDs. 

In this paper we explore the possibility that inner rims are 
much higher when the gas and dust temperatures are computed sep- 
arately and self-consistently. A detailed study of inner and outer 
thermo-chemical disc structure s together w i th em ission lines in 
Herbig Ae discs is presented in Ka mp et all d2010h . In this paper 
we focus on the inner disc structure of well-mixed discs and the 
images and spectral signature of inner rims in the infrared. 

This paper is organized as follow. We briefly describe the 
ProDiMo code and recent additions to the code such as the SED 
and image calculation in Sect. [2] We proceed with a presentation of 
our fiducial disc model, the disc parameters that are varied, and the 
fixed-structure disc model in Sect. [3] We then present and discuss 
the results of our simulations when the gas and dust temperatures 
are equal (thermal-coupled models) and when they are computed 
self-consistently (thermal-decoupled models), varying several disc 
parameters in Sect. [4] Finally, we conclude about the inner disc 
structures and their effects on the SEDs in Sect. [5] 



2 PRODIMO CODE DESCRIPTION 

ProDiMo is designed to compute self-consistently the (1 + 1)D 
hydrostatic disc structure in thermal balance and kinetic chemi- 
cal equilibrium. The disc can be active or passive depending on 
the choice of the viscosity parameter a. We first provide here a 
brief description of ProDiM o. Interested readers are referred to 
Woitke, Kamp & Thi (2009) to find detailed explanations of the 
physics implemented in the code. We continue by giving the ratio- 
nal for the "soft-inner edge" surface density profile implemented 
in ProDiMo. We finish by explaining the computation of Spec- 



tral Energy Distributions and images, which is a new feature of the 
code. 

2.1 General description 

The code iterates the computation of the disc density, gas and dust 
temperature, and chemical abundance structure until successive it- 
erations show less than 1% change. 

The iteration starts with a given hydrostatic disc structure 
computed from the previous iteration or from an assumed struc- 
ture. The dust radiative transfer module computes the dust temper- 
atures and the disc local mean specific intensities J„(r, z). The spe- 
cific intensities are used to calculate the photochemical rates. The 
2D dust radiative-transfer module of ProDiMo has been bench- 
marke d against other co des that use other methods such as Monte- 
Carlo dPinte et al. 2009). The continuum radiative transfer includes 
absorption, isotropic scattering, and thermal emission. The grains 
are assumed spherical. The lack of anisotropic scattering prevents 
us to calculate completely accurate images. In case of anisotropic 
scattering, one expects the far upper part of the disc rim to be dim- 
mer (backward scattering) than the near lower part (forward scatter- 
ing). However, isotropic scattering is sufficiently accurate for gen- 
erating precise SEDs. 

The grains have sizes between a m i n and a max and follow 
a power-law size distribution with index p. The absorption and 
scattering efficiencies are computed using Mie theory for compact 
spherical grains. Grains of different sizes are assumed to have the 
same temperature. This limitation in our model does not prevent the 
comparison between self-consistently computed gas temperature 
models (thermal-decoupled models) and models with T gas = Td ust 
(thermo-coupled models). In this paper, the dust composition and 
size distribution are constant through out the disc. We a dopt the in- 
terstellar dust optical constants of Laor & Draine J 1993b for amor- 
phous grains since we are only interested in modelling featureless 
broad SEDs. 

Once the dust temperature Td ust (r, z) and the continuum 
mean specific intensities J v {r,z) are known, the gas temperature 
and chemical concentrations are consistently calculated assuming 
thermal balance and kinetic chemical equilibrium. The chemistry 
network includes 71 gas and ice species. The chemical reactions 
comprise photo-reactions (photoionization and photodissociation) 
with rates computed using cross-sections and specific intensities J„ 
calculated by the 2D radiative transfer module, gas phase reactions, 
gas freeze-out and evaporation (thermal desorption, photodesorp- 
tion, and cosmic-ray induced desorption), and H2 formation on 
grain surfaces. The chemical rates are drawn from the UMIST 
2007 database dWoodall et alll2007l) augmented by rates from the 
NIST chemical kinetic website. The photo-cross-sections are taken 
from ryan Dishoe ck et al. I d2006h . The H2 formation on grain sur- 
faces follows the prescription of ICazaux & Tielensl d2004l) with the 
most recent grain surface parameters. 

The gas can be heated by line absorption and cooled by line 
emission. The line radiative transfer is computed via the escape 
probability formalism. The gas is also heated by interactions with 
hot photoelectrons from dust grains and Polycyclic Aromatic Hy- 
drocarbons and with cosmic rays. Thermal accommodation on 
grain surfaces can either heat or cool the gas depend ing on the tem - 
perature difference between the gas and the dust (Tielens f2005l) . 
Thermal accommodation dominates at high densities and optical 
depths, driving the gas and dust towards the same temperatures. 

Finally, the gas temperature and molecular weight determine 
the local pressure at each grid point, which is used to modify the 
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Figure 1. The upper right and left panels show the dust and gas temperature structures respectively for the fiducial model. The dust temperature contours are 
overplotted for T(j ust = 50, 100, 500, and 1000 K in both panels. The location of Ay = 5 and 50 are also shown. The entire outer disc at R > 100 AU is 
located at Ay < 5. The lower-right panel show the disc density structure for the entire disc (up to Rout = 300 AU). The gas temperature over dust temperature 
ratio for the fiducial model in the lower-left panel. The gas can be more than ~ 5 times warmer than the dust in the upper disc atmosphere but T gas = Td us t 
at Ay > (1 — 5). The ratio T g!lB /T^ uat is explicitly truncated at 5 in the figure. 



disc structure according to the vertical hydrostatic equilibrium. The 
new disc structure is compared to the previous one to check for 
convergence. 

2.2 Soft Inner edge 

The radial dependence of the vertical column density is as- 
sumed to deviate from a power-law in the innermost layers, 
and is implemented according to the "soft-edge"-description of 



IWoitke, Kamp & Thill20cM (see Section 3.1). We assume that the 
gas has been pushed inwards by the radial pressure gradients 
around the inner edge, and has spun up according to angular mo- 
mentum conservation, until the increased centrifugal force balances 
the radial pressure gradient. The procedure results in a surface layer 
in radial hydrostatic equilibrium at constant specific angular mo- 
mentum where the column density increases gradually from virtu- 
ally zero to the desired value at the point where we start to apply the 
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Table 1. Disc parameters. When a parameter has multiple entries, the values in bold correspond to the values of the fiducial model. 
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2.2 Mq 
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effective temperature 


en 


8600 K 


disc mass 


M disc 


10~ 2 , 10" 3 , 10" 4 Mq 


disc inner radius 


Rin 


0.5, 1 , 10 AU 


disc outer radius 


Rout 


300 AU 


vertical Column density power law index 


e 


1, 1.5,2 


gas to dust mass ratio 


S 


100 


dust grain material mass density 


Pdust 


2.5 g cm -3 


minimum dust particle size 


^min 


0.05 Atm 


maximum dust particle size 


^max 


10, 50, 200 Atm 


dust size distribution power law 


V 


3.5 


H2 cosmic ray ionization rate 
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1.7 x 10- 17 s- 1 


ISM UV field w.r.t. Draine field 


X 


0.1 


abundance of PAHs relative to ISM 


/pah 


0.01, 0.1 


a viscosity parameter 


a 


0.0 
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Figure 2. Inner rim density structure and photospheric height (ylv=l con- 
tour) for the fiducial model ("soft-edge" and T gas computed by thermal 
balance). 



power-law. The thickness of this "soft" edge results to be typically 
a few 10% of the inner radius Ri n . 



2.3 Spectral Energy Distribution calculation 

Based on t he results of the co ntinuum dust radiative transfer so- 
lution fsee l Woitke et alj d2009h . Sect. 4) we have developed a new 
ProDiMo module for SED and image calculation based on formal 
solutions of the dust continuum radiative transfer equation along a 
bundle of parallel rays that cross the disc at given inclination angle 
i with respect to the disc rotation axis. The rays start at an image 
plane put safely outside of the disc. The ray direction ft and the 
3D origin of the image plane p& are given by 



(1) 




.£? i(r 




. _ _ a_ = 200 um 

a m „ = 50 urn 

a„„ = 10 urn 



io~'L 
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X [|lm] 
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1000.0 



Figure 3. The dust opacity computed by Mie theory for the three values 
of Umax (10, 50, 200 /im) and a power-law dust grain size distribution 
(index 3.5). A line showing an opacity with a A -1 opacity law is added for 
comparison. 



and 

pfj = R~ft, 



(2) 



where R is a large enough distance to be sure that the plane cannot 
intersect the disc. The plane is aligned along the two perpendicular 
unit vectors 



(3) 




and 




(4) 



The rays are organized in log-equidistant concentric rings in the 
image plane, using polar image coordinates (r,9). The 3D starting 
point of one ray xfy, with image coordinates x = r sin(#) and y = 
r cos(#), is given by 
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xf,^p^ + xnt + yn^ (5) 

From these ray starting points xq, the radiative transfer equation 
is solved backwards along direction —~n, using an error-controlled 
ray integration scheme as described in Woitke, Kamp, & Thi (2009, 
Sect 4.2). The source function and opacity structure in the disc have 
been saved from the previous run of ProDiMo's main dust radia- 
tive transfer. The resulting intensities at the location of the image 
plane, I v (r, 9), are used to directly obtain monochromatic images. 

To retrieve the spectral flux under inclination i, we integrate 
over the image plane as 

Fv = h J l »( r > e ) rdrde © 

where d is the object's distance. We use Ne = 160 and N r = 200 
rays. Due to the log-spacing of the concentric rings, the method 
safely resolves even tiny structures originating from the inner rim. 



3 DISC PARAMETERS 
3.1 Fiducial model 

We modelled circumstellar discs around one typical Herbig Ae star 
with effective temperature of 8600 K, mass of 2.2 M©, and lumi- 
nosity of 32 L©. In this study, we did not try to fit the SED of 
any particular object. The input stellar spectrum is taken from th e 
PHOENIX database of stellar spectra dBrott & Hauschildtll2005h . 

The discs are illuminated from all sides by a low interstellar 
medium UV flux (\ = 0.1) in addition to the stellar flux. We chose 
a low flux because previous studies do not take interstellar UV flux 
into account. The stellar flux dominates over the standard interstel- 
lar flux in the inner disc by orders of magnitude. 

The abundance of PAH (/pah) in discs controls the gas tem- 
perature and the disc flaring, which in turn affects the dust tem- 
perature. The exact value of /pah can only be constrained by a 
simultaneous fit to the broad SED and the PAH features for a spe- 
cific object. We chose an arbitrary low abundance (1%) and tested 
the effects on the SED for the fiducial model with /pah=0. 1. 

The only parameters that were allowed to change are the disc 
total mass (M d i sc = 10~ 4 , 10" 3 , and 10~ 2 M ), the surface 
power-law index (E = E D x (r/r )~ e with e — 1.0, 1.5, and 2.0), 
the disc inner radius (7?i n = 0.5, 1 , 10 AU), and the maximum 
grain radius (a max = 10, 50, 200 /jm). 

We chose e = 2 as our fiducial value for the surface density 
profile index. Surface density profiles as steep as r~ 2 for the in- 
ner d isc are a signature of self-gravitating discs dRice & Ar mitage 
2009). However, we consider discs with maximum disc mass over 
star mass ratio of 4.5 x 10~ 3 . A recent derivation of the Sol ar Neb- 
ula al so gives a steep surface density profile inside 50 AU jDeschl 
120071) . Shallower decreases in surface density (e = 1.5) are more 
commonly observed for the outer disc (R > 50 AU). 

The density falls off at R[ n following the "soft-edge" model 
described earlier. The minimum grain size a m in is kept at 0.05 /im. 
The power-law index of the dust grain size distribution follows the 
standard interstellar value of -3.5. A power-law index of -3.5 re- 
sults from grain-grain collision theory dDohnanvil 1969i) , We mimic 
grain growth by varying the maximum grain radius a max from the 
fiducial value of 50 to 200 /im. We also studied the effect of smaller 
grains by setting the maximum grain size to 10 (jm. 

The choice of the inner disc radius at 0.5 AU was dictated 
by disc accretion and dust grain physics. A disc is truncated at the 



co-rotation radius in the absence of a magnetic field or of a low- 
mass companion, either star, brown dwarf, or giant planet, orbiting 
close to the star. The co-rotation radius r Iot = GAf*/i> 2 ot is de- 
fined as the distance from the central star where the Keplerian disc 
rotates at the sam e speed than the star photosphere v Iot = v, (e.g., 
IShu et alj fl994). Assuming a photospheric rotation speed v t = 
50 — 150 km s _1 , typical of young stars ( I van den Ancker et al.l 
1 19981) . we obtain r rot ~ 0.08-0.76 AU for M, = 2.2 M Q . Dust 
grains do not exist above their sublimation temperature, which 
is around 1500 K for silicate grains. The dust sublimation ra- 
dius is the distance from the star beyond which dust exists and 
can be estimated by r d = ^/Qr(L, + L acc )/167rcr/r 2 ub where 
Qr = Qabs(a, r*)/Q a b s (a, T sub ) the ratio of the dust absorption 
efficiency at stellar temperature T, to its emission efficiency at the 
dust sublimation temperature T su b, a is the mean grain radius and 
a is the Stefan constant. For large silicate grains (a ^ 1 pirn), Qr is 
relatively insensitive to the stellar effective temperature and close to 
unity because most stellar radiation lies at wavelengths shorter than 
the grain r adius. For smaller grains, the val ue of Qr is significantly 
increased dMonnier & Millan-Gabetll2002h . For one micron grains, 
the sublimation radius in our fiducial model is rd = 0.19 AU as- 
suming L acc << L* (passive disc) and Qr ~ 1 (large grains). 
For grains as small as 0.05 in radius, we obtain Qr ~ 20 and 
rd = 0.88 AU. The inner disc location at 0.5 AU in our mod- 
els is within the theoretical range of inner disc radii. Moreover 
grains as small as 0.05/im can survive at 0.5 AU. Large luminosity 
(L» > 100 Lq) fast rotators (usini > 200 km s" 1 ) like 51 Oph 
(with Z/»= 26OI5Q L© and u*=267±5 km s _1 ) may have a in- 
ner dust-free gas-rich reg ion between the co-rotation radius a nd the 
dust sublimation radius jThi et alj|2005l ; iTatulli et alj|2008h . An- 
other example is AB Aur where the corotation radius is smaller 
than the sublimation radius for the small grains but not for the large 
grains. Observational evidence of an inner gas-rich dust-po or re- 
gion has been found for this object jTannirkulam et al1l2008l) . 

We also chose to simulate discs with inner radius at 1 and 
10 AU to model typical gaps possibly created by planets. The outer 
radius R ou t was set at 300 AU. By default passive discs are mod- 
elled (i.e., a = 0). Since we focus on young discs the standard 
interstellar medium value of 100 for the gas-to-dust mass ratio was 
adopted for all models. The gas and dust were assumed well-mixed 
with no dust settling. The common model parameters are summa- 
rized in TableQ]with the parameters for the fiducial model printed in 
boldface . For each parameter set, the code was run in the thermal- 
coupled mode (Td us t = T gas ) and thermal-decoupled mode in a 
50 x 50 non-reg ular grid. An exam ple of the location of the grid 
points is given bv lKampetai] f2010). All SEDs were computed for 
a source located at the typical distance of 140 pc at two inclinations: 
face-on (0 degree) and 45 degree. 

The stellar properties of our fiducial model differs from earlier 
studies (T off = 10000 K, M»=2.5 M , and L«=47 Lq). We plot in 
Fig.|9]the density structure, photospheric height, near-infrared im- 
age, and SEDs for two inclinations using the earlier studies prop- 
erties. The disc structures and SEDs are different for the two sets 
of stellar parameters. The hotter disc atmosphere around the 47 L© 
star flares much more than the disc around the lower luminosity 
star. The infrared luminosity is increased but the shape of the SED 
does not differ much as the stellar luminosity is increased. 



3.2 Fixed-structure model 
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Figure 4. Spectral energy distribution computed for the various models. The upper left panel shows the fiducial model (Af ( jisc=10 — 2 Mq) with decreasing 
surface power index e and in dot-dash purple line the SEDs computed for the fixed disc structure models with the flaring index p=1.25 and e=2. The upper-right 
panel shows the SED for the fixed-structure model with e = 2 and three values of the flaring index p. The upper-left panel is for discs seen face-on. All the 
other panels are for discs seen at 45 degree. The other panels show the effect of varying the disc total mass (10 — 2 , 10 — 3 , and 10 — 4 Mq), the grain upper size 
limit (a max = 10 fim and a max = 200 fim), the inner disc radius (0.5, 1, and 10 AU), and the effect of adopting a "sharp-edge" density profile compared 
to a "soft-edge" density profile. SEDs in solid lines are computed in the case T gas is computed by thermal balance while SEDs in dash lines are computed 
assuming T gas = T^ust- The SED of the fiducial model is shown in red. The input stellar spectrum is plotted in light blue. The spectral sampling is much 
higher for the stellar spectrum than for the SEDs. The flux in the UV and optical is the sum of the direct stellar flux and the stellar flux scattered by the disc. 
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Figure 5. Inner rim density structure for the fiducial model but with different disc structure prescriptions as indicated above each panel. 



Table 2. Disc parameters for the fixed-structure model. 

reference scale height Ho 0.72 AU 
reference radius Ro 10 AU 

flaring index p 1.2,1.25,1.3 



We compare the results of the hydrostatic disc models to the results 
of a fixed- structure model. The disc height H is parametrised by 
the functional H(R) — Ho(R/ Ro) p , where Hq is the reference 
scale height at reference radius Ro and p is the flaring index. The 
scale height is used to compute the density nn(R, z) at height z by 
nu(R,z) — no(R)exp( — (z/H) 2 ). The value of the parameters 
are given in Tableland were chosen such that the SED matches 
that produced by the fiducial model with "soft-edge" and T gas com- 
puted by thermal balance. We vary the value of the flaring index p 
between 1.2 and 1.3 to study its effects on the SEDs. All the other 
parameters are the same than for the fiducial model. 



4 RESULTS AND DISCUSSION 

In this section, we discuss the effects on the disc gas density struc- 
ture and SED when we vary a few parameters in the case of equal 
gas and dust temperature and when they are computed indepen- 
dently. We first show the structure and SED of the fiducial disc 
in Sect. 14.11 The SEDs and images of the fixed-structure disc is 
discussed in Sect. 14.21 We continue by addressing the effects on 
the disc structure and SED when we use the "soft-edge" versus 
the "sharp-edge" inner rim prescription (Sect. l4~3l >. when we vary 
the surface density profile (Sect. I4.5t . the disc mass (Sect. 14.6b . 
the maximum grain size (Sect. I4.7K and the inner disc radius 
(Sect.|48]l. 



4.1 Fiducial model 

Figure [T] contains some results for the fiducial model. The dust 
and gas temperature distribution are shown in the upper two pan- 
els while the lower-right panel is the disc density structure. The 
disc areas where the gas and dust temperature are decoupled are 
emphasized in the lower-left panel. Figure [3] shows the dust opac- 
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Figure 6. Image of the inner disc at 1.98 fjm for the fiducial model (Afdisc=10 Mq, e = 2) for the two possible combinations of inner edge prescriptions 
and two ways to compute the gas temperature and for the fixed-structure disc structure. In the fixed structure case, the scale-height is 0.72 AU at radius 10 AU 
and the flaring index is 1.25. The disc is seen at an inclination of 45 degree. The inner radius was set at 0.5 AU. We can see a second rim at 1-2 AU, which is 
the most pronounced in the case of soft-edge and gas temperature computed by thermal balance. The emission from the second bump contributes significantly 
to the total flux. The fixed-disc model has a flaring index of p = 1.25. 



ity (scattering and absorption) for three maximum grain sizes (10, 
50, and 200 /im). The opacity is dominated by scattering at wave- 
lengths below 10 fim and by absorption at longer wavelengths. 

In the thermal-coupled models, the dust temperature struc- 
ture is divided into a hot inner rim, a warm upper layer, and a 
vertical isothermal interior. The lower-left panel in Fig. Q] shows 
the ratio between the gas and dust temperature. The gas and 
dust temperature are equal at vertical optical depths Ay > 1 for 
a 10~ 2 Mq disc. The 2D dust continuum radiative-transfer re- 
sul ts confirm the analytical tw o-zone decomposition popularized 
bv lChiang&Goldreichl Jl997h . The SED of the fiducial model is 
shown in solid-red lines in Figure|4] Model images of the inner disc 
at 1.98 fim and 4.96 fim are shown in Fig.[6]and Fig.|7]respectively. 
The flux at 1.98 fim comes mostly from the inner rounded rim. The 
region just behind the rim is deprived of stellar photons, is cooler, 
and thus does not emit strongly (shadowing effect). 

The gas temperature remains above 1000 K in the disc atmo- 
sphere up to a few AU. Together with the decrease in the vertical 
gravitational pull, they explain the emergence of a secondary 
density bump at ~ 1.3 AU. The presence of the hot "finger" in 
disc atmospheres is typical of disc mo dels that compute the gas 



and d ust temperature independently tercolano. Clarke & Drake] 
120091 : iGlasseold. Meiierink & Naiital 120091 ; IWoitke. Kamp & Thil 



20091: iNomura & Millail 120051: Ikamp & pullemonol |2004 
Jonkhei d. Faas. va n Zadelho ff & van Dishoeck 2004). 

The shape of the disc photospheric height (Ay=l con- 
tour) is shown in Figure [2] The height at 0.6 AU is located 
at z/r ~ 1, The curvature exhibits the same shape than in 

models where grains of different sizes sublime at different radii 

i ir— — i 

(Isella& Natta 2005) or with density-dependent sublimation tem- 
perature jTannirkulam et alj|2007t) . The curvature of the inner rim 
is caused by gas angular momentum conservation and pressure gra- 
dient and is independent on the dust physics. Likewise the location 
of the inner rim is set by hydrodynamic constraints (disc truncation 
or presence of a companion) and not by dust sublimation physics 
(see Sect. l3.lt . However we have assumed that the gas and the dust 
are well-mixed. Therefore, the gas and dust grain density structures 
may be independent if grain-growth and settling occur. 



4.2 Fixed-structure model 

The inner disc structure for the Ho = 0.72 AU and flaring index 
p =1.25 fixed-structure model is shown in the upper-right panel of 
Fig. [5] The inner rim height is ~ 0. 1 AU. The SEDs that correspond 
to the fixed-structure models are drawn for inclination (p=1.25) 
and 45 degree (p=1.2, 1.25, 1.3) on the two top panels of Fig. [4] The 
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Figure 7. Same as Fig.[6]but the images were generated at 4.96 (im 



flaring index p impacts on the amount of warm dust grains in the 
disc atmosphere that emit beyond 100 /im. The SED of the fixed- 
structure model with flaring index p =1.25 matches relatively well 
the SED of the fiducial model. The images at 1.98 and 4.96 /im for 
the p = 1 .25 model are shown on the upper right panel of Fig.|6]and 
Fig. 13 

When viewed exactly face-on, the SEDs at 3 /mi of a fixed- 
stru cture disc show very weak near-IR fluxes, which arise from the 
rim jMeiier et al.ll2008f) . The cause of this feature is that the pro- 
jected inner rim emission area is null when the disc is seen face-on. 
In contrast, the hot-dust emitting area in the soft-edge models with 
rounded inner rim show much less dependency on the viewing an- 
gle. 

At 45 degree inclination the flux at 3 (im matches the fluxes 
from all the other disc-structure models. The flux comes mostly 
from the edge facing the star, which cannot be seen at null incli- 
nation. The SED between 8 and 100 /im is extremely sensitive to 
the value of the flaring index p because the more flare a disc is, the 
more stellar photons are intercepted by the outer disc, which pre- 
dominately emits in the far infrared. The short wavelength flux (3-8 
/im) is mostly sensitive to the height of the rim, which is constant, 
and changes only a little when we vary p. 

SEDs alone cannot be used to differentiate between the in- 
ner disc structures. On the other hand, images at short wavelengths 
(see Fig.[6]l show differences between the models. The emission at 
1.98 /im of the fiducial model is concentrated in the rounded in- 
ner rim with contribution from the second rim. The emission of the 



fixed-disc model is spread over the first AU. Discriminating disc 
models by fitting simultaneously the SED and images has already 
been use d successfully to study the disc around the T Tauri star 
IM Lupi dPinte et alj|2008l) . 



4.3 Testing the effect of "soft" versus "sharp" edge 

Current disc models assume that the disc density decreases mono- 
tonically from the inner edge to the outer radius ("sharp-edge" 
model). We study the effect of adopting our "soft-edge" model 
compared to "sharp-edge" model on the disc structure and SEDs. 
Figure [5] shows the density structure for the fiducial disc parame- 
ters assuming T gas = Td ust on the left panel and of T gas computed 
by gas thermal on the middle panel. Those two density structures 
contrast with the two upper panels in Fig. [8] 

First the rim is extremely narrow in the "sharp-edge" models. 
On the other hand, the main effect of the 'soft-edge' model is to 
create a rounded rim surface whatever the surface density profile 
(see Fig. [5}. The density gradients at the rim are clearly seen in 
the density structure plots. The projected ro unded rim appears as 
an ellipse, similar to other theoretical studies Jlsella & Nattall2005l : 
iKama et alj|2009l) . 

In the "sharp-edge" models, the strength of the 3 /im excess 
emission depends on the inclination of the disc with respect to the 
observer. This dependency stems from the thinness of the rim and 
was already found in previous studies of thermal-coupled models 
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Figure 8. Inner rim density structure for the fiducial model (JWdisc=10 Mq) with surface power index e =1.5 and e =1. The left panels are the models 
with Tgas = 7d ust ■ The right panels are the models with T gas computed by gas thermal balance. These models show a much taller rim and also a secondary 
rim at ~ 1.2 AU. 



jMeiier et alj|2008b . On the other hand, the 3 /im excess depends 
much less on the inclination for the "soft-edge" models (see upper- 
left panel of Fig.|4j. 

Secondly, there is no secondary "bump" in the structure at 
~1.2 AU, even when T gas is computed by thermal balance. The 
presence of the second bump in the "soft-edge" decoupled model 
can explain the difference in fluxes between the models in the 4 and 
40 micron region in the lower-right panel of Figure[4] The optically- 
thick rim in the "sharp-edge" models prevents a large amount of 
photons to reach the disc surface at radius 1—1.5 AU (the so-called 
"shadow"), therefore suppressing the rise of the second bump. The 
bump exists because the gas stays at temperatures > 1000 K up to 
4 AU at the surface while the dust temperatures are only > 300 K. 



The bump creates a large emitting area for the warm dust, which 
translates into strong emission in the 4^10 fim region in the SEDs 
(See Fig.|7]and upper-left panel of Fig.|4j- 

In the "sharp-edge" models, the mid-infrared flux is depressed 
compared to the "soft-edge" models because of the lower disc 
height behind the rim. 

4.4 Effects of PAH abundance 

The P AHs are the main heating of the agents of the gas 1 Ka mp et al.l 
120101) . A disc with / P ah=0.1 is warmer than a disc /pah=0.01. A 
warmer disc is more extended vertically and intercepts more radia- 
tion from the star (Fig.[l2]upper-left compared to Fig.[8]upper-right 
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Figure 9. Inner disc density structure (upper-left panel), photospheric height (upper-right panel), near-infrared image (lower-left panel), and comparison 
between the SED for the fiducial model (lower-right panel) but with M„=2.5 Mq, T c ff =10,000 K, and L*=47 Lq. We show the L* = 32 Lq stellar spectrum 
only. The red-dashed line in the upper-left panel encloses the region where the gas temperature is higher than 1000 K. Notice the large amount of scattered 
light in the blue part of the visible range in the SEDs caused by the extremely flaring inner disc atmosphere for the more luminous star. 



panel) although the shape of the photospheric height does not vary 
significantly (Fig. Q2] upper-right panel). The overall emission in 
the continuum and in the gas lines is higher (Fig. 1121 lower-right 
panel). We also show an image at 1.98 in the lower-left panel 
of Fig. [T2] 



4.5 Effects of varying the surface density profile 

We show the results of varying the surface density power-law in- 
dex from e = 1,0 to e = 2.0. All other parameters are kept con- 
stant. The density structure for the thermal-coupled (Td ust = T gaa ) 
and thermal-decoupled (T gas from thermal balance) models are dis- 
played in Fig.[8]for the inner disc on the right and left respectively. 



The Spectral Energy Distribution for the models are shown in the 
upper-left panel of Fig. [4] 

The structure of thermal-decoupled disc models differ sub- 
stantially from their thermal-coupled counterparts. The gas is 
heated by gaseous and dust grain photoprocesses that convert ultra- 
violet photons into fast-moving electrons, which in turn share the 
energy to the gas (mostly atomic and molecular hydrogen). In the 
upper disc layers, the density is too low for efficient gas-grain ther- 
mal accommodation and the gas and dust are thermally decoupled. 
The gas mostly cools by line emissions, which become quickly op- 
tically thick, while dust grains cool by optically thinner continuum 
emission. As a result, the gas remains at higher temperature than 
the dust grains. 
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Figure 10. Inner rim density structure for the A/ t ji sc =10 3 Mq (upper panels) and M^ SC =10 4 Mq (lower panels). The left panels are the models with 
Tgas = 7d ust . The right panels are the models with Tgas computed by gas thermal balance. 



Below the disc atmosphere, the ultraviolet flux is attenuated 
and thermal accommodation between gas and dust grows. The gas 
there is mostly molecular and frozen onto grain surfaces. Finally, in 
the highly extinct midplane thermal accommodation dominates and 
drives towards equal gas and dust temperatures. The dust tempera- 
ture structure is not affected by the change in density structure as 
there are still two layers: a warm upper layer and a vertical isother- 
mal interior. 

The inner rim is much higher when the gas temperature is de- 
termined by detailed heating and cooling balance. The SED reflects 
relatively well the amount of hot grains and thermal-decoupled 
models show a larger amount because of higher rims. 

In the thermal-coupled model, the rim is much low er and does 
not bl ock photons to reach the area behind the rim. lAcke et all 
d2009h studied the SED from Herbig Ae discs using a Monte-Carlo 



code coupled with hydrostatic equilibrium assuming equal temper- 
atures. They found good matches to observational data when they 
artificially rise the height of the rim, consistent with our results 
when the gas thermal balance is computed. 

The second bump at 1-2 AU that appears in the thermal- 
decoupled models are not present in the r gas = Td us t models 
(the left panels in Fig. [8). As already mentioned in the discus- 
sion of "sharp-edge" versus "soft-edge" models, the second bump 
is caused by the gas that is hotter than the dust at the disc surfaces. 
A rounded optically thinner rim helps more scattered light to reach 
the disc surfaces compared to a narrow optically thicc rim. 

The existence of a puffed-up rim is a necessary but not suffi- 
cient condition for strong 3 /J,m excess. The surface area of the rim 
has to be large, i.e. the rim has to be high enough. In the seminal 
paper on disc inner rim by iDullemond et all j200 lh . the height of 
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Figure 11. Inner rim density structure for the a m ax = 10 /im and a max = 200 (im. The left panels are the models with T ga s = 7d us t. The right panels are 
the models with T gas computed by gas thermal balance. 



the rim to the rim distance from the star ratios were found to be of 
the order of 0.1-0.25. We find that the rim in the case of e = 2.0 
is 0.25 AU high at 0.5 AU. Here we have defined the disc height as 
the location of the disc where gas density is 10 6 cm~ 3 . A higher 
rim is need ed in our modelling to pr oduce strong emission at 
because the lDullemond et all ( 1200 lh model assumed a single dust 
temperature at 1500 K for the rim while in our models only the dust 
in the rim at the midplane reaches 1500 K. In the upper-left panel 
of Fig. Q] we can see a vertical decreasing dust temperature gradi- 
ent from the midplane (1500K) to the surface (< 1000 K) at the rim 
radius. 

The mid-infrared continuum emission at a specific wavelength 
reflects the amount of dust grains at a temperature roughly given 
by Wien's law. As found in previous studies, discs flare in the outer 



disc in the thermal-coupled models. The rim deprives the disc at- 
mosphere behind the rim of UV photons. This shadowing effect 
has been invok ed to explain SEP of disc w i th low 30-100 /xm 
emission flux jPullemond & Dominild |2004 iMeiier et all l2009t 
IVerhoeffetaT]|2010h . This shadowing effect is less pronounced if 
dust scattering is taken into account. Dust scattering ensures that 
some photons reach the "shadowed" area behind the rim. 

The main effect of warmer upper layers in the decoupled mod- 
els is puffed-up and strongly flaring upper layers from the inner rim 
to the outer radius. Although the inner rim is twice higher than in 
the thermal-coupled case, the extra dust opacity in the rim is not 
sufficient to attenuate the UV flux responsible for the gas heat- 
ing behind the rim. Denser gas are found higher in the disc and 
more grains are raised to higher temperature, resulting in slightly 
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Figure 12. Inner disc density structure (upper-left panel), photospheric height (upper-right panel), near-infrared image (lower-left panel), and comparison 
between the SED with /pah=0 01 and /pah=01 for the fiducial model (lower-right panel). The red-dashed line in the upper-left panel encloses the region 
where the gas temperature is higher than 1000 K. 



stronger 10-50 /im continuum flux. In general our hydrostatic disc 
models flare strongly. 

4.6 Effects of varying the disc mass 

In this series of models, the surface density profile is kept at 2.0 
and we vary the disc mass (Mdi sc =10 -4 , 10 -3 , and 10 -2 Mg). 
The structure and SEDs for this series are shown in Fig. 1101 andBl 
respectively. The near-IR flux is weakly affected by the disc mass. 
On the other hand, the flux in the 10-50 /im region is sensitive to 
the total dust mass and thus decreases with decreasing total disc 
mass (the gas-to-dust mass ratio is kept constant at 100), consistent 
with the findings of lAcke et ail d2009h . From the SEDs in Fig. 0] 
the disc mass is the most important parameter controlling the near- 



to far-IR flux ratio. As the opacity decreases with wavelength (see 
Fig. O, the flux depends less and less on the disc geometry and 
more and more on the total dust mass, which scales with the gas 
mass if a constant gas-to-dust mass ratio is assumed. 

4.7 Effects of varying the maximum grain size 

The maximum grain size has a relatively weak influence on the 
inner disc structure (Fig.ll U. Therefore the flux at 3-5 /im does not 
change significantly with increasing dust grain upper size limit (see 
Fig. 0. On the other hand, the effects on the 10-70 fim shape of 
the SED are significant. The bigger the grains are, the cooler they 
are and the emission is weaker as testified by the higher opacity 
per unit mass for the small grains [3] In addition, big grains have 
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Figure 13. Inner rim density structure for the Ri n =l (upper panels) and 10 AU (lower panels). The left panels are the models with T gas 
panels are the models with T gas computed by gas thermal balance. 
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less surface area per volume and therefore the photoelectric heating 
efficiency decreases. The disc is cooler and flares less, which results 
in less intercepted stellar photons. The combined effects concur to 
decrease the 10-70 flux when grains are big. Beyond 100 /im, the 
grain emissivity depends strongly on the presence of larger grains. 



At wavelength larger than their sizes, grains are efficient emit- 
ters. Therefore the flux drop at wavelengths longer than 70 /im is 
less pronounced for big grains. In summary, increasing the max- 
imum grain radius results in weaker emission at 30-70 fim but 
stronger emission at longer wavelengths. 



4.8 Effects of varying the inner disc radius 



Fi gure |4l and 1 1 3 1 show the SEDs and inner disc structure for Ri n = 1 
and 10 AU. A round rim structure is present for all inner disc radii. 
Since the hottest grains are cooler, the peak of the near-infrared 
bump has shifted to ~ 4 fj,m in the Ri n = 1 AU model. The spec- 
tral signature of the inner rim has disappeared for the Ri n = 10 AU 
model. The presence of the puffed-up rim does not depend on the 
actual value of the inner radius but only rims close to the star have 
hot enough dust to emit strongly in the near-infrared. 
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5 CONCLUSIONS 

We modelled the structure, Spectral Energy Distribution, and in- 
frared images of protoplanetary discs around Herbig Ae star and 
analysed the difference between the assumption of equal gas and 
dust temperatures and when both temperatures are computed inde- 
pendently. We also compared the hydrostatic disc model SEDs to 
the SEDs of a disc with prescribed density structure. 

The disc structure is governed by the gas pressure support, 
which in turn depends on the gas molecular weight and temper- 
ature. The height of the inner rim in models with calculated gas 
temperature exceeds those where the gas and dust temperature are 
equal by a factor 2 to 3. Higher rims result in large emitting areas 
and thus slight larger near-infrared excess. Our treatment of the in- 
ner disc density fall-off ("soft-edge") results in rounded inner rims 
consistent with near-infrared interferometric studies. The discs also 
show a second density bump that manifests itself as stronger emis- 
sion between 3 and 30 fim. The flux beyond 30 pim is mostly sensi- 
tive to the disc mass. The maximum grain radius affects weakly the 
SEDs. The effect of gaps is to remove hot dust grains which emit 
predominately in the near-IR. As a result, the flux in the near-IR is 
suppressed. The shape of the SED from 3 to 100 /im cannot be used 
to discriminate between the inner disc structure models. Together 
with the SED, images with high spatial resolution may be used to 
differentiate between disc models. 

Our study stresses the importance of understanding the inter- 
play between the gas and the dust in protoplanetary discs. This in- 
terplay shapes the disc structures, which in turn control the shape 
of the SEDs. 
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